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We study a generalization of a Pais-Uhlenbeck oscillator for fermionic variables. Next, we consider 
an ensemble of these oscillators and we identify a particular case of the Myers-Pospelov model 
which is relevant for effective theories of quantum gravity. Finally, by taking the advantage of this 
connection, we analyze, for this model, the unitarity at one loop order in the low energy regime where 
no ghost states can be created on-shell. This energy regime is the relevant one when we consider 
the Myers-Pospelov model as a true effective theory coming from new space-time structure. 



I. INTRODUCTION 

In recent years, there has been growing interest in 
quantum field theories containing higher time deriva- 
tives. In part because they arise naturally as effective 
corrections from candidate fundamental theories such as 
string theory [JJ, quantum gravity [5] and noncommu- 
tative geometry [3J. Moreover, attempts to incorporate 
such terms have also been given in the context of dark 
energy [4], Lorentz invariance violation [Ml], radiative 
corrections [SI H] perturbative iteration [TUj and regular- 
ization [TT] . 

Many years ago, a quantum mechanical system with 
higher time derivatives was proposed by Pais and Uhlen- 
beck (P-U) 12J. This system provides a generalization of 
the harmonic oscillator by including a fourth time deriva- 
tive. The first problem that one faces is the unbounded 
character of the spectrum |13j : sometimes, however, con- 
strained systems can turn around the problem, see [14] 
and there has been also other proposals [15] . There are 
some redefinitions which avoid these problems but in turn 
they introduce states of negative norm. Lee and Wick 
(L-W) generalize these ideas to relativistic quantum field 
theories [TB]. These undefined metric theories have the 
advantage of being finite in the sense that these higher 
time derivatives regularize the divergences. The prob- 
lem of unitarity, which comes out with undefined metric, 
could be controlled by a prescription due to Cutkosky 
[TT]. Recently, these models have been used as an ex- 
tension of the standard model of particles which fixes 
the hierarchy problem |18j . However, models incorporat- 
ing operators of mass dimension 5, which violate Lorentz 
symmetry, might also contain higher time derivatives. 
These theories due to Myers and Pospelov (M-P) [5] are 
supposed to arise from quantum gravity effects at high 
energies. These two last examples, are part of the in- 
creasing interest on higher time derivatives in quantum 
field theory and, therefore, it is compelling to focus on 
the technical and physical issues concerning these kind 
of theories. 

In this work, we propose to construct a field theory 
as ensembles of Pais-Uhlenbeck oscillators in momentum 



space, instead of harmonic oscillators of the standard 
relativistic case. We find that this new theory is close 
related to the Myers-Pospelov fermionic sector. Taking 
advantage of this relation, we explore the unitarity of the 
model by modifying the Cutkosky prescription. 

All the theories with higher time derivatives seem to 
have an unavoidable problem because they content neg- 
ative norm states, at least if there is no gauge symmetry 
that prevents these modes from appearing in the spec- 
trum 14^ The fermonic sector of the Myers-Pospelov 
model also has this feature, and these ghost modes ap- 
pear even at tree level in perturbative theory, when the 
initial state have a total energy greater than the mass of 
one of these modes. However, this model has to be seen as 
an effective model, only valid for energies below certain 
scale, beyond which new physics is described through 
a more fundamental theory. Hence, lack of unitarity 
and probabilistic interpretation problems coming from 
the ghost states are irrelevant because the description 
of physics beyond that threshold with this model is not 
justified. Nevertheless, the model must be well defined 
below the threshold, and it must be capable of giving 
sensible predictions at this low energy regime. But, uni- 
tarity could be spoiled out due to virtual processes even 
though no on-shell ghosts can be created at low enough 
energies. This problem was addressed by Cutkosky et 
qL[T7] in the context of relativistic Lee- Wick theories, as 
mentioned above. They used a convenient prescription 
to define the propagators in order to keep unitarity or- 
der by order in the loop expansion. We will follow the 
same philosophy to define the propagators in such a way 
that the Myers-Pospelov model preserves unitarity below 
the threshold scale. It is worth noting that in the work 
by Cutkosky the Lorentz invariance is and important in- 
gredient to be considered. However, in our present case, 
relativistic symmetry only plays a secondary role which 
will serve to fix part of the prescriptions which connects 
to the standard relativistic theory at very low energies 
compared with the threshold. 

The layout of this work is the following. In the second 
section, we review the Pais-Uhlenbeck oscillator and we 
generalize it to fermions. In the third section, we con- 



struct an ensemble of P-U fermionic oscillator and we 
identify it as the fermion sector of M-P model. In sec- 
tion four, we analyze the unitarity of the Myers-Pospelov 
model at one loop order, in the low energy regime by us- 
ing the structure of P-U oscillators. Finally, we give the 
conclusions and final comments. 



II. THE ISSUE OF HIGHER TIME 
DERIVATIVES 

In this section, we introduce the P-U model and we 
generalize it to the fermionic case which may be less fa- 
miliar, and hence it will be studied with more details. 



where a+, a + , a_ , a_ are the standard creation and anni- 
hilation operators. The second term produces arbitrary 
negative energy states as can be seen by acting a_ on the 
empty wave function (defined by a+^o = a-^o = 0) 



$o = N exp 
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An alternative proposal to quantize would be to redefine 
the vacuum a+<&' = a/_& = 0. Then the states ip n = 

a_ wo have positive energy, but are ghost states when n 
is odd. Hence, this procedure stabilizes the theory but it 
could spoil unitarity when interactions are introduced. 



A. The Pais-Uhlenbeck model 

The P-U oscillator is, basically, the standard harmonic 
oscillator with an additional higher time derivative term. 
To be more precise its equation of motion is 



fj<l 



(4) 



tA = o, 



(i) 



where q^' is a fourth order time derivative and g can be 
considered a small coupling constant. The equation of 
motion is obtained from the Lagrangian: 
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(2) 



This system can be seen as two standard harmonic oscil- 
lators by means of the change of variables, 

q + = {8 2 t - k 2 _)q, 

q_ = (df k\)q. (3) 

The Lagrangian with the new variables is 



jpu 



I 2 2 
k, Q, 



\1- 



-klqt, 



(4) 



with k± = 4(1^ yj\ — 4gw 2 ) positive frequencies de- 
pending on g and u>. This Lagrangian corresponds to two 
standard harmonic oscillators with one of them having a 
relative minus sign respect to the other. This last fact 
makes a great difference with the simple sum of two os- 
cillators when we add interactions or quantize the theory 
Classically, the equations of motion are just the same as 
two standard harmonic oscillators, however when adding 
interactions the system becomes unstable due to the un- 
boundedness of the energy. This problem persists under 
the usual quantization. However, by a redefinition of the 
vacuum state, the quantum mechanical problem become 
stable. This last procedure leads to unavoidable negative 
norm states. These ghosts states could render the the- 
ory non unitary when interaction are considered. To be 
more precise, following the canonical formalism we write 
the Hamiltonian as: 



Hpjj — k + a + a+ 



fc_a!a_ 



-(fc+-fc_), (5) 



B. Fermionic extension of the P-U model 

Now we consider the fermionic version of the previous 
model with equation of motion: 



g%l> + iip — ujip = 0. 



(7) 



As before, this is the fermionic harmonic oscillator (g = 
0) plus a higher derivative term. Without loss of gener- 
ality we can assume that uj and g are positive constants. 
This equation of motion comes from the Lagrangian 



Lp = gt/jip + v^ty ~ ojipij). 
Let us consider the new variables 

V»_ = a(id t -u + )ip, 



ip_ = a(—idt — <*>+)ij>, 



and 
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where a = 
are defined by 



and the positive frequencies u>± 
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where ui- > u + > 0. The Lagrangian in terms of this 
new variables is 



L = i/)Uid t - u+)ip+ - i>-(idt - uj-)ip- 



(12) 



which can be reduced to the original Lagrangian pi), by 
replacing Eqs. ([9| and (10 1 on it. 

Now, the usual variation of the above Lagrangian gives 
the equations of motion 



(id t 
(id t 



w+)^+ = o, 

LJ-)lp- =0. 



(13) 



As in the bosonic case the system is decoupled into two 
standard oscillators, but with a global minus sign for the 



-0_. Following the canonical quantization for fermions, 
one can check that the canonical conjugate momentum 
for ip + is iip + , but for ip_ is —iip_. Hence, this gives the 
anticommutators in the new variables 



where we have made use of the anticommutation relations 



{tl>+,il>\.} = 1, 



(14) 



Then, the operators ^>_, ip_, ip+ and ip\ are standard 
creation and annihilation operators, respectively, but ■0_ 
and ip_ create and annihilate states of negative norm. 

The original variables can be written in terms of these 
operators as 



V> 



ip- - VM 



{l-Aguy/v 



ijj = —i 
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so the Hamiltonian turns out to be 

Hp = w+i/j + i/j+ — u)-ip_i/j- 



(15) 



(16) 



This expression can be compared to the one for the P- 
U Hamiltonian. They are the same expression, except 
for an irrelevant constant. However, this system corre- 
sponds to a four state system, and hence, the spectrum 
is bounded below, and the vacuum is normalizable. In 
fact, the true vacuum is the state annihilated by ip_ and 
;/'+, namely, ^A-'&o — V'+'&o — where the explicit form 
of the wave function for the vacuum has been calculated 
in the appendix. 

Moreover, the rest of the states are $i = ip + $o, with 
energy u+, $2 = 4>-®o with energy w_ and $3 — 
'0-V'+^ > o with energy lo + + cj_. Nevertheless, the $2 and 
$3 operator correspond to ghost states. In the appendix 
we follow the Schroedinger quantization procedure and 
the interpretation is the same as above. 

Therefore, in the fermionic case we just have one alter- 
native to quantize the theory, in contrast with the bosonic 
case where we can either have unstable spectrum or in- 
definite norms in agreement with the analysis with Lee 
and Wick [16]. 

Finally, we calculate the temporal ordered Green's 
functions which will be useful in later perturbative con- 
siderations. To do this, let us consider the usual defini- 
tion 

G F ab {t) = 0(*)<O|i&«(tWj(O)|O) - 9(-t)(0\4(0)%(t)\0) 
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(17) 



where a,b 
we have 



±. By solving the equations of motion (13) 



i>a(t) 



'MO). 
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Replacing this in the Green function expression and tak- 
ing into account that ^ a (0)|0) =0we arrive at 



G F ab {t) 



ae 



l e{t)5 ab 



(19) 



(14). This Green function satisfies the equation 



(id t - uj a )G F b (t) = ai5(t)5 ab . 



(20) 



We can use an integral representation with the following 
prescription 



G F b {t) = aiS ab 



du> 



27r uj — io a + ie ' 



(21) 



where the prescription e > is used. 

Note that we are dealing with a (0+1)— dimension 
model where there is no particles and antiparticles. In 
a more realistic (3 + 1)— dimension case, particles and 
antiparticles are present, and the usual prescription for e 
is positive for particles and negative for antiparticles. 



III. FERMIONIC PAIS-UHLENBECK 

ENSEMBLE AND THE MYERS-POSPELOV 

MODEL 

In this section we propose a generalization of the 
fermionic P-U model to a quantum field theory. This 
generalization may be performed in different ways. How- 
ever, we will use the following approach: we consider 
the standard field theory as an ensemble of decoupled 
harmonic oscillators in the momentum space, with and 
additional term containing higher time derivatives, as we 
did in the previous section. 



A. Covariant free field theory as ensembles of 
harmonic oscillators 

To fix ideas, let us review how we can go from the 
standard bosonic and fermionic harmonic oscillator, to a 
scalar and Dirac free field theories, respectively. 

First, we consider an ensemble of bosonic harmonic 
oscillators whose dynamical variables are labeled by a 
vector p with frequency up, so that the Lagrangian of 
this ensemble is a sum of the Lagrangian of the individual 
oscillators over all possible values of p. In the continuum 
limit the sum becomes an integral, and we can write the 
Lagrangian as 



L 



1 



d^p^+uj}) 



(22) 



where 4>p is the dynamical variable of the oscillator. In 
order to make contact with a relativistic theory we as- 
After 



yfp 2 + m 2 . After a Fourier transform 
we end up with the free scalar field theory 



sume that cup 



L 



d 3 xct) S (D + m 2 )(j) S , 



(23) 



where cj)g(t) is the Fourier transform of 4>p(t). 

The fermionic case is less direct. In order to recover 
the free Dirac field, we will need four fermionic oscillators 



at each point in the p vector space. We will label this 
fermion variables as ipf~ where s will be the spinor index 
and i will specify particle or antiparticle modes. Then, 
the Lagrangian at each point p is 



L p = E ^ idt 



'GW&. 



(24) 



Again, the ensemble of this oscillators will be given 
through the integral over all possible values of p. To get 
the relativistic theory, we assume that the frequencies 
have the form 



u!(p) = ti^p, 



(25) 



with e,; — ±1 for particles and antiparticles, respectively, 
and we write the new variables 



Mt) 



i 



]T<(p)V^(t), 



(26) 



where w\{p) = u s (p) is the amplitude of the free wave 
function that is solution of the Dirac equation for parti- 
cles of spin s and momentum p and tu|(p) = V ' S (~P) the 
free wave function for antiparticles of spin s and momen- 
tum —p. In terms of these basis and making use of the 
orthogonality relation 



> 



wV(?)w s j (p)=2u; i! 8 rs 8V, 



and the completeness relations 



, w i (f) w f (f) = Wp + €i h D (p) , 



(27) 



(28) 



where ho(p) = a ■ p + to/3 and a — 7°7 and /3 = 7 are 
the 4x4 Hermitian Dirac matrices, the Lagrangian turns 
out to be 



-L/p ■ 



^(idt 



h D {p)^p- 



(29) 



Hence, the Lagrangian of the ensemble, which is the in- 
tegral of the above expression, can be written in terms 
of the Fourier transform variables ips and "ijjg = ip s (3 as 



L = 



d 3 p 
(2tt) 3 f ~ 



d 3 xipz(id ■ 7 — ra)i\)z. (30) 



Summarizing, this procedure shows us that given an 
ensemble of oscillators, with frequencies depending on 
free parameters, we can find a space that we call ir-space 
where the ensemble is described as a relativistic standard 
free field theory. 



B. An ensemble of P-U harmonic oscillators: 
Fermionic Myers-Pospelov model 

Now, we come to the main part of the section. As it 
was mentioned before, we will consider an ensemble of 
P-U harmonic oscillators and see what theory produces. 



We set out from Eq. ( 24 ) and ( 25 ) , but we add a higher 
time derivative term 



L p = E ^dt + idt ~ e^)<- 



(31) 



Then we consider the variable i/y as defined in Eq. ( 26 ) 
and performing a Fourier transform the Lagrangian has 
the form 



L 



d 3 p 

(2tt)3 p 



d 3 Xipg(id ■ 7 + gd\ - m)ip s (32) 



This model, arising from an ensemble of P-U harmonic 
oscillators, can be identified with fermionic sector of the 
well known M-P model, which is an effective Lorentz in- 
variance violating field theory, containing quantum grav- 
ity effects. To clarify this claim, we introduce a vector 
n = (1, 0, 0, 0), and the former Lagrangian can be written 
as 



L = d xips(i0 + gi/i{n ■ d) — rnjips- 



(33) 



This is the Myers and Pospelov fermionic sector with 
.92 = and an isotropic background with 5-dimension 
operators. This model has been widely studied in the 
context of quantum gravity effects O [19] . 



IV. UNITARITY 

Indefinite metric theories have potential physical prob- 
lems which may spoil out any predictive ability and even 
worse, they would lack of a sensible physical interpreta- 
tion [50] . The source of these problems is the failure of 
the unitarity of the time evolution operator. Many years 
ago, Cutkosky proposed a prescription to keep the uni- 
tarity at perturbative level for nonhermitian but Lorentz 
symmetric theories [16 . This procedure was successfully 
performed for the so called Lee- Wick models. 

In this section, we attempt to fit similar prescription 
in order to prove the unitarity for the M-P fermionic 
sector at one loop perturbative level, and for energies less 
than the ultraviolet scale of the theory. This modification 
must take into account the Lorentz invariance violation 
intrinsic in the M-P model. As it was mentioned in the 
introduction, we are only interested in this low energy 
regime because it is the relevant one for a theory which 
must be seen as an effective theory parametrizing new 
features of space-time structure which may come from a 
more fundamental theory. 

The first thing to do is to write down the relevant prop- 
agators. In the M-P theory the propagators, in terms of 
fundamental fields are complicated, and it is more conve- 
nient to express them in terms of the P-U modes. These 
are the P-U modes 
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(34) 
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FIG. 1. We show the diagrams that contain ghosts and there- 
fore those that can spoil unitarity at one loop order. In the 
first figure we show the diagram corresponding to the ampli- 
tude M^ , in the second, the amplitude M_ + , and in the last 

figure the amplitude corresponding to M 



in terms of the Fourier transform fields 
ijj(t,p) = J d 3 xip(t,x)e lp ' x . Now these new variables cor- 
respond to four P-U oscillators with frequency uj{p) = 
CiUif. As we have seen in section II this oscillators are 
equivalent to two fermionic oscillators with frequencies 



u r ±M 



2g ' 



for each P-U oscillators labeled with i and r. The ± signs 
refers to the signs of the norm for those modes. Now, it 
is straightforward to calculate the propagators, 

<?;$(!,% = (0|T(^ a (i,x)^(0)|0), (36) 

corresponding to these modes by looking at the expres- 
sion at the end of the section II. Namely the Fourier 
transform of this propagators are, 



G Fab( UJ 'P)= a 



V -L , 

2u> p - u - uliij?) + id? 



(37) 



where 8f will correspond to some suitable prescription 
to fit later. In order to recover the standard Lorentz 
symmetric limit we must impose the prescription 8^ = 



—82 = e > 0. To see this, the standard Feynman propa- 
gator can be related with the P-U propagators by 



A F (u,p) = J2 G F+>>P)"t° 



(38) 



Then, in the limit g — > we have uj s +i (p) — > eiCJp, and 



making use of the completeness relations ( 28 ) we obtain 
the standard relativistic expression 



A F (p) = 



2wp 



Wp + h D (p) 



h D {p) 
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jr 



(39) 



Now, we will prove that with a suitable prescription for 
8~ we can keep the unitarity of the theory at one loop 
perturbation, at least. Let us write the evolution opera- 
tor in terms of the scattering matrix, 

U = 1 + iS. (40) 

Then, the unitarity implies, 

UU ] = l + i{S-S^) + SS^ = 1. (41) 

Namely, we must prove that — i(S — 5>) = SSr. 
Perturbatively, we must show that 

2ImM(in|out) = ^ 7W(in|Phys)X(Phys|out),(42) 

Phys 

where M. are the perturbative elements of the scattering 
matrix and the sum is over all possible physical states 
allowed by the conservation of energy and momentum. 

We are interested in showing this equality for in and 
out states for particles with energies less than the ul- 
traviolet scale j-. At tree level, due to the energy and 
momentum conservation, it is impossible to create ghosts 
states and then, the contribution is zero at both sides of 
the previous equation, at first order in perturbation cor- 



rection. The first chance to break down the equation (42 ) 



(35) is at one loop order; however, for the same reason above 



mentioned, the right hand side of equation ( 42 ) is zero 



Then, we must check that the imaginary part of M is 
zero at the first radiative correction. 

It is sufficient to analyze the integrals corresponding to 
the diagrams in the Fig. 1, which have problematic ghost 
internal lines. To simplify, we will work in the center of 
mass frame, with total energy Q and total momentum 
P = 0. The integrals contributing to this graphs are, 



M ab (fl) 



^2 / dwd 3 p 



i,i,i',j' 



tr (g£"(« - n/2,f)Gn/ J '(" + 0/2, *0) (43) 

where a is the coupling constant which introduces an 
interaction and a, b = ±. 



By using the orthogonality property in the Eq. ( 27 ) the 



sum of the integrands of these diagrams are proportional 
to the functions, 

r a6 M = V — p^ 7^— x 

^(jj- S2/2 - Uai(P) + 10" 



1 



(44) 



w + Q/2 - u; bi (p) + iSl 

The pole structure of these integrands is shown in Fig. 2 
and Fig. 3, and their residues at the poles are 

1 



Res r Qb = 
Res r ab = 

U(,i+fi/2 


U>ai 


— Ubi — 
1 


(45) 
(46) 
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- Wai + 



Any prescription to avoid the poles is equivalent to a path 
integration in the w-plane C. 

For cop < j-, we choose the prescription 6± = — #J = 
5 > 0. When we integrate u over the real axis, this 
prescription leave the poles 1~ below and 2~ above. This 
is equivalent to the contour C shown in Fig. 2. The 
conjugate of this integration is given by the contour C* . 
Hence, the imaginary part of the integral is proportional 
to the closed contour integration over C — C* , which is 
proportional to the sum of the residues. 

Thus, for hxiM we have, 

f T duo ex Y^ e^Res V (w), 

JC-C* „.,„w_<±§ 



-»-*M 



E« 



1 



1 



= 0. 



And for Im(M.-\ + M |_), we have 

f (r+_ + T + _)duj ex V e l Res (r+_ + T_+) 
Jc-c* i , w±i±n/2 



E- 



uj + i — uj—i — Q uj-i — ui+i + Q 



0, 



where e^ stands for signs ±1 of the closed contours ori- 
entation that surround the poles. 

For internal momentum above the critical scale, top > 
i , some poles become complex, as it can be seen in Fig. 
3. Now the integral over the real axis of to can be closed 
around the lower half-plane, picking up the residues at 
the poles in this half-plane. The complex conjugate inte- 
gration will be closed above the upper-half plane, picking 
up now the poles in this half-plane. Thus, the imaginary 
part of the integral with complex poles is given by a con- 
tour enclosing these poles. This prescription is equivalent 
to take the poles corresponding to 1 + moving up and the 
ones corresponding to l - moving down. Then, it can 
be seen that the prescription given above corresponds to 
take the integration path C and C* shown in Fig. 3, which 
are smooth deformations of those paths in Fig. 2. 

Now, the closed contour integrals take the residues in 
the same way as before, and their sum vanishes again. 
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FIG. 2. We show the pole structure for uip < l/4p, the no- 
tation i stands for particles and antiparticles with positive 
or negative norms (±). The crosses correspond to the shifted 
poles in the integrals of the amplitudes at one loop order. 
The contour C corresponds to the prescription described in 
the text, C* to its conjugate. The vertical line corresponds 
to the critical region above which the poles become complex. 
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FIG. 3. The pole structure above the critical region ujp > 
l/4g. The poles corresponding to 1 + and 1~ collapse to each 
other and they become complex. According to our prescrip- 
tion the pole 1 + go up and 1~ go down in the critical vertical 
line. The integration contours C and C* are smooth defor- 
mation of those of the Fig. 2. 



Thus, we have found a particular prescription to main- 
tain the unitarity at one loop order. 

We have used a somehow arbitrary interaction, how- 
ever, the procedure followed above can work with other 
general interactions. 

Finally, let us note that there is some ambiguity in 
choosing the prescription which keep the unitarity. This 
freedom could be used to maintain the unitarity at higher 
radiative corrections. Nevertheless, this is far from being 
clear in this analysis. 



V. FINAL COMMENTS 

Summarizing, in this work we have found a connection 
between an emsemble of fermionic P-U oscillators, and 
the well known M-P model in the fermionic sector. This 
fact has allowed us to prove the unitarity of the M-P 
fermonic sector at one loop order in perturbation the- 
ory at low energies where the theory is applicable, even 
though, it contains higher time derivative terms. This 
statement is the most important result in this study, be- 
cause it permits to have sensible predictions and physical 
interpretations. 

Our study shows that for Lorentz invariant violat- 
ing theories with higher time derivative terms - with 
fermions - it is possible to maintain the unitarity at first 
order in radiative corrections. 



The prove of the unitarity at higher radiative correc- 
tions, is left as a future work. Also, it would be interest- 
ing to do phcnomenological predictions from this theory. 
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Appendix A: Schroedinger Quantization of the P-U 
Fermionic Oscillator 



In this appendix we show an alternative procedure to 
quantize the P-U fermionic oscillator. This formalism 
gives us an explicit form of the vacuum state, and shows 
how the negative norm states appear. 

The following Hamiltonian 



H = —gipip + unpip, 



and the non vanishing anticommutators 



(Al) 



9 9 9 



reproduce the fermionic P-U equations of motion pre- 
sented in section II. 

Using the Schroedinger representation in terms of 



Grassman variables and their derivatives, we have, 

f i 9 i - 

gdip 2g 

; i d i , , . „. 

tP=-w7-^- ( A 3) 

And the Hamiltonian in this representation is, 



H : 



+ uj 1 
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(A4) 



The first and second lines of this equation commute to 
each other, so it is very easy to find the eigenf unctions, 



*o =ip, 

(l-4g«)4 

$ 2 = - ^ , e- 

(l-4gu)4 

*3 = i>, 



s/l-igu. 
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- tprp 



Ei =w+, 

E 2 =u-, 

E 3 = lu + + u- , 



(A5) 



where w± are defined in section II. For positive <?, the 
state $o corresponds to the lowest energy state, and then 
it represents the vacuum of the theory. Now, it is clear, 
that this is a four state system with bounded energies. 

Regarding the normalization, we can define the Berezin 
measure as J dip dip ipip = +1, the scalar product of wave 
functions F(ip, ip) and G(ip, ip) is 

(G\F) = f dip dip GF. 

Then we can see that 

(fco|$o> = <*i|*i> = l. (A6) 

However, 

($ 2 |$ 2 ) = ($ 3 |$ 3 ) = _i (A7) 

are states with negative norm. This agrees with the re- 
sults in section II. 
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